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Abstract 

We present an analytical treatment of the dissipative-stochastic dy- 
namics of a charged classical particle confined bi-harmonically in a plane 
with a uniform static magnetic field directed perpendicular to the plane. 
The stochastic dynamics gives a steady state in the long-time limit. We 
have examined the orbital magnetic effect of introducing a parametrized 
deviation (77 — 1) from the second fluctuation-dissipation (1I-FD) relation 
that connects the driving noise and the frictional memory kernel in the 
standard Langevin dynamics. The main result obtained here is that the 
moving charged particle generates a finite orbital magnetic moment in 
the steady state, and that the moment shows a crossover from para- to 
dia-magnetic sign as the parameter r\ is varied. It is zero for 77 = 1 that 
makes the steady state correspond to equilibrium, as it should. The mag- 
nitude of the orbital magnetic moment turns out to be a non-monotonic 
function of the applied magnetic field, tending to zero in the limit of an 
infinitely large as well as an infinitesimally small magnetic field. These 
results are discussed in the context of the classic Bohr- van Leeuwen theo- 
rem on the absence of classical orbital diamagnetism. Possible realization 
is also briefly discussed. 

PACS numbers: 05.40.-a, 05.10.Gg, 75.20.-g. 

The Bohr-van Leeuwen (BvL) theorem [T]-[3] on the stated absence of or- 
bital diamagnetism for a classical system of charged particles in equilibrium has 
been one of the surprises of physics [3] . The static external magnetic field ex- 
erts a Lorentz force on the moving charged particle, acting at right angle to its 
instantaneous velocity (v). While such a gyroscopic force does no work on the 
particle, it does induce an orbital cyclotron motion that subtends an amperean 
current loop. The magnetic field associated with this current loop is expected 
to be non-zero, and directed oppositely to the externally applied magnetic field 
- the Lenz' law. Hence the expectation of a finite orbital diamagnetic moment 
q/2c(r x v) [5]. Yet, as is known well, the partition function for a classical sys- 
tem in equilibrium turns out to be independent of the applied magnetic field, 
thus giving a zero orbital magnetic moment. And this has been the surprise 
[4]. (The field independence of the classical partition function derives simply 
from the fact that the classical partition function involves integration of the 
canonical momentum over an infinite range for any given value of the conjugate 
coordinate, and thus the magnetic vector potential (A) entering the Hamilto- 
nian through minimal coupling (p — > p — |A) gets eliminated through a trivial 



1 



shift of the momentum variable. This shift is, however, not allowed for a quan- 
tum system because of the canonical non-commutation involved there. Hence 
the stated quantum origin of the orbital magnetic moment in equilibrium - the 
Landau orbital diamagnetism 6J. A remarkably heuristic real-space explana- 
tion for the vanishing of the classical orbital moment was first suggested by 
Bohr pQ in terms of a cancellation of the orbital diamagnetic moment of the 
completed amperean orbits (Maxwell cycles) in the bulk interior by the para- 
magnetic moment subtended by the incompleted orbits skipping the bounding 
surface of the system in the opposite sense. This 'edge current 1 has a large 
arm-length, or leverage and, therefore, can effectively cancel out the bulk dia- 
magnetic moment. The cancellation has indeed been demonstrated graphically 
for a simple planar geometry [7] . This real space-time picture is consistent with 
the zero orbital magnetic moment following from an exact analytical solution 
of the classical Langevin dynamics with a white noise describing the motion of 
the charged particle confined harmonically in two dimensions, with a uniform 
static magnetic field applied perpendicular to the plane [8]. Here, the steady 
state (t — >• oo) orbital moment indeed vanishes for the given potential confine- 
ment (owing to the spring constant k of the harmonic potential, providing a soft 
boundary). Interestingly, this null result persists in the limit k — » 0, provided 
it is taken after taking the limit t — > oo. This suggests that in this case the 
stochastic particle dynamics has had time enough to sense (i.e., be influenced 
by) the confinement (k). On the other hand, the moment survives to a non-zero 
value if the order of the two limits is interchanged. (The effect of these so-called 
Darwinian limiting processes is also manifest in the case of the quantum version 
of the above Langevin treatment [9 )■ It is to be noted, however, that in the 
case of quantum Langevin equation the orbital moment tends to zero as the 
Planck constant is formally reduced to zero, i.e., in the classical limit, for which 
the noise term reduces to a classical white noise which is consistent with a local 
Stokes friction constant - in accord with the second fluctuation-dissipation (II- 
FD) theorem (TU1 [TT]. This reasonably suggests to us that it may well be the 
constraint of the second fluctuation-dissipation relation that forces the orbital 
magnetic moment to vanish in the classical case. This is further supported by 
our recent numerical simulation [T2]. Motivated by these observations, we have 
carried out an exact analytical calculation of the orbital magnetic moment of 
a charged particle confined bi-harmonically in two dimensions with a uniform 
static magnetic field applied normal to the xy-plane, but now with the proviso 
that the stochastic driving force (noise) is a sum of two uncorrelated noise terms 
- an exponentially correlated term and a delta-correlated term - and there is a 
parametrized (jj — 1) deviation from the II-FD relation. Interestingly now, we 
do obtain a non-zero orbital magnetic moment in the infinite time limit - in 
the steady state. Moreover, the sign of the orbital magnetic moment turns out 
to show a dia- to para-magnetic crossover as the parameter rj is tuned through 
rj = 1, where the moment vanishes. In the following we present an exact analyt- 
ical treatment of this dissipative-stochastic system and discuss the results that 
follow. 

Consider the classical dissipative-stochastic dynamics of a particle, of charge 
— e and mass m, which is confined bi-harmonically in the ir?/-plane in the pres- 
ence of a uniform static magnetic field B applied perpendicular to the plane. 
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The governing stochastic (Langevin) equations are 



mx{t) = -kx(t) - J ^I e - (t -*' )/tc + T S(t - Ht'W - ^tf(*) + fx(*) 



where £f(t) is the noise term with (&(t)) = 0, and 

(£*(W)) = Sijmk B T ^l e -l*-*'l/*. + 2r, 7o( 5(i - t')) 



(lb) 



(2) 



and we are interested in the long-time limit t — ► oo. Here i = x,y, and the 
angular brackets (• • •) denote average over realizations of the two un-correlated 
noise terms - one a delta-correlated (white) noise and the other an exponentially 
correlated noise with a correlation time t c . This sum of a white noise and 
an exponentially correlated noise, we believe, is the simplest non-Markovian 
gaussian process allowed by Doob's theorem [15] , Further, (jj — 1) parametrizes 
deviation from the II-FD relation as noted above. 

It is convenient to introduce here the quantities f2o = y/k/m (harmonic os- 
cillator circular frequency) , — = J7 C (cyclotron circular frequency) , and ^ = 7 
(the frictional relaxation frequency). We further define the following dimension- 
less parameters r = yt (dimensionless time), and ujq = ^°-,w c = %*■ (dimension- 
less circular frequencies). 

Following now the 'Landau trick', the two coupled Langevin equations for 
the real displacements x(r) and j/(r) as functions of the dimensionless time r, 
can be conveniently combined into a single Langevin equation for the complex 
displacement z(r) = x(t) + iy(r), giving: 

z( T ) = -u%z(t) - jT (L e -{r-r')/r c + J^ S ^ T _ T 'fj i( T ') dr ' + iu c z{ T ) + g{r) 

1 ' (3) 

with (g(r)) = and (g(r)g*(r')) = (i e -(— + { 2 V )^S(r - r')) . 

Note the complex conjugation (*) that we have introduced in (g(r)g* (t 1 )) 
above, as the same will be needed in subsequent calculations. Also, we have 
changed over to the dimensionless time parameter (r) , but have retained the 
same symbols for the dynamical variables without the risk of confusion. Inas- 
much as the particle motion along the uniform magnetic field normal to the 
xy-plane decouples from that in the xy-plane, the present model equally well 
describes a 3-dimensional system. The orbital magnetic moment can now be 
re- written as 

(M(r)) = glm<z(r)i*(r)). (4) 
It is convenient now to introduce the Laplace transform 

/>oo 

z(s) = / e~ ST z{r)d,T, and z(s) = sz(s), (5) 



with the initial conditions z(Q) = = z(0) (without loss of generality, as we are 
interested in the steady state (long-time limit r — > oo). We obtain straightfor- 
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wardly 

( s3 +(i + ?-- C ) 52 +(-o 2 + i + ^-^) S +^) 



m = 7 7 9(s)(l + r c s)/r c 



and 



|( a ) = g(l + r eS )g( S )/r e 

(* 3 +(^ + ?-- C )^+(-o 2 + i + ^-^) 5 + ^)' 

where we have used i(s) = sz(s), with the initial conditions z(0) = z(0) = as 
noted above. In order to inverse-Laplace transform the expressions above, it is 
convenient to introduce the factorized denominator 

D(s) = (s^ + (^ + ^-i. c )s^ + U + ^ + ^- l ^-) S+ ^ 

V \ T c 7 J \ T c "/T c T c J T c 

= (s - si)(s - s 2 )(s - s 3 ), (8) 

where Sj(i = 1,2,3) are the three roots of the cubic denominator D(s). These 
roots are readily obtained following the Cardano procedure. We can then write 
z(s) and z(s) as partial fractions 



^) = Et^-^m, (9) 

and 

^) = Etzt^' (10) 

S^ 

where {ai} are given as solutions of the set of equations 
ai + a 2 + a 3 = 

ai(s2 + S3) + 02(53 + si) + 03 (si + s 2 ) = -1 

1 , x 

ais 2 s 3 + a 2 s 3 s 1 + a 3 sis 2 = — , (11) 

and similarly, {A{\ are given by 

A1+A2 + A 3 = 1 

Ai(s 2 + s 3 ) + A 2 {s 3 + si) + A 3 (si + s 2 ) = 

^is 2 s 3 + A2S3S1 + A 3 s lS2 = 0, (12) 
In terms of the above, we now have the inverse Laplace transforms as 

z(t) =Va t f e s ^-^/^g(T r )dT' (13a) 
i Jo 

z{t) = V / T e s ^- T "^^g(r")dT" . (13b) 



and 
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With this, the orbital magnetic moment in Eq. (4) turns out to be 



mc/ \ 7r ( 



i,j=l,2,3 



8i(r-T')/T e Sj(t-t")/t c —\t'-t"\/t c 



+ 2, / 2^£ e ^(r-r')/r CeS -(r-r")/r^ (r / _ ^ 
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(14) 



where Im denotes the imaginary part. Straightforward integration gives the 
t — y oo (steady-state) limit for the orbital magnetic moment as 



M <~> = (£) (w) Im S a,A 

V ' c 7 8,7=1,2,3 



1 



I I „ 70 
rr f - 2?7— r, 



7 V Si + Si 



• (15) 



After some simplification, the above expression for the orbital magnetic mo- 
ment M(oo) reduces to 



M(oo) 

( efcgT \ 



=- 2(?? -i)^im x: 



1 
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i,j=l,2,3 



S l 



(16) 



In Fig. 1, we have plotted the limiting steady-state value of the dimensionless 
orbital magnetic moment fj, against the dimensionless applied magnetic field 
j3(= ^f-), where rj parametrizes the II-FD violation and is varied over the 
range 0.5-2.0. As we are interested here mainly in the matter of principles, we 
have made a simple choice for the dimensionless parameters involved, namely 
^- = uiq = 1 and ^ = 0.5 (the strength of the white noise relative to the 
exponentially correlated noise). 

As we observe the field-induced orbital magnetic moment is clearly non-zero 
in the steady state. There is a crossover from the paramagnetic to the diamag- 
netic sign as the parameter rj is tuned from rj = 0.5 to rj — 2.0. This crossover 
is a surprise. The magnetic moment is zero for rj = 1, which corresponds to 
the canonical II-FD consistent (equilibrium) state. Hence no violation of the 
BvL theorem. The orbital magnetic moment is obviously zero for zero magnetic 
field; but, not so obviously it tends to zero for large magnetic field as well. The 
latter behaviour can, however, be understood from the following, namely that 
the radius/frequency of the cyclotron orbit tends to zero/infinity as the applied 
magnetic field is made infinitely large 0]. The fact that the orbital magnetic 
moment can be paramagnetic in certain range of the II-FD deviation parame- 
ter rj is significant in that, unlike diamagnetism, it leads to a positive feedback 
for a collection of charged particles in such a classical system - it can give an 
enhancement of the orbital paramagnetism. 

Physical realization of such a classical system in the laboratory is admittedly 
somewhat demanding. One needs to create a dilute (highly non-degenerate) gas 
of charged particles (e.g., electrons/holes) at sufficiently high temperatures, and 
confined on a mesoscopic scale in the presence of a static magnetic field. The 
temperature has to be high enough so as to wash out the quantum effects, 
namely the discreteness of the quantized level spacings owing to the mesoscopic 
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Figure 1: Plot of dimensionless magnetic moment /i against dimensionless mag- 
netic field /3 for four different values of r\ that parametrizes deviation from the 
II-FD relation. It clearly shows a dia- to para-magnetic crossover as r\ is varied 
through r) = 1. Also, the moment tends to zero in the limit of zero as well as 
large magnetic field /?. 

confinement. Now, the II-FD violating parameter (77 — 1) necessarily requires a 
non-equilibrium steady-state condition. This is the real problem for an experi- 
mental realization. One is tempted to think that such a non-equilibrium steady 
state may be induced through a noisy laser excitation, e.g., the Kubo- Anderson 
non-Markovian noise [141 115] , of charged particles confined in an optical tweezer 
(16j . There is, however, a problem here involving the energy injected by the 
laser, its dissipation in the system and the associated rise of temperature. In 
fact, one necessarily needs to have a two-temperature configuration. Thus, e.g., 
a possible physical realization of our model can, in principle, have charged parti- 
cles in contact with a thermal reservoir at one temperature, while another type 
of neutral particles is in contact with the reservoir at a different temperature. 
Then particle collisions will ensure a stationary (steady-state) non-equlibrium 
state with flow of energy between the reservoirs. This then according to our 
model calculation should give a non-zero orbital magnetic moment in an exter- 
nally applied magnetic field |17j . 

Given that classically a static magnetic field does no work on a moving 
charged particle, our model calculation giving a non-equilibrium steady-state 
solution in the presence of a static magnetic field could lead to some insightful 
molecular dynamical (MD) simulations when appropriately thermostatted [18] . 
Also, it has a significant bearing on the work related to generalized fluctuation- 
dissipation theorem for steady-state systems j!9]. 

In conclusion, we have presented an exact analytical treatment of a classical 
dissipative-stochastic model system in a uniform static magnetic field, which is 
found to give a finite orbital magnetic moment in the steady state. Interestingly, 
we find that there is a crossover from the diamagnetic to the paramagnetic sign of 
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the magnetic moment as function of a parametrized deviation from the second 
fluctuation-dissipation relation. We think that these results do complement, 
rather than violate the classic Bohr-van Leeuwen theorem. 

The author would like to thank A. A. Dcshpandc and K. V. Kumar for many 
discussions. The author would also like to thank the referee for constructive 
criticisms. 
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